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1 Introduction 



Integrable quantum spin chains are quantum mechanics models which can be solved exactly. Cer- 
tainly, among them, the most studied is the spin | Heisenberg spin chain (also called the XXX 
chain) [1] solved by H. Bethe [2] using his celebrated ansatz. However, in recent applications as con- 
densed matter experiments [3-5] or string theory (for recent reviews, see [6,7] and references therein), 
more involved models show up. These models deal with higher rank Lie algebras such as gl(N) (the 
XXX chain being associated to gl(2)) and/or higher spin (i.e. higher dimensional representation) on 
each site, see e.g. [8-18]. They also include spin chains with impurities [19-21] and alternating spin 
chains [22-25]. In most papers, each model is investigated on its own. Obviously, in opposition to 
these case by case studies, it would be of very much interest to develop a general framework which 
allows one to solve all these models simultaneously. It is provided by the study of an 'algebraic 
version' of spin chains which encompasses a generic spin chain where each sites may carry differ- 
ent representations of gl(N). The resolution of this general model started with the computation of 
the Bethe equations in [12,26,27]. Let us note that a very similar approach has been developped 
in [28, 29] and [30] where an 'algebraic version' of Bethe vectors is computed, in term of L-operator 
(first case) or Drinfeld currents (second case). Then, the vacuum state of this model was constructed 
in [31] in the thermodynamical limit (i.e. when the length L of the chain tends to infinity) when the 
representations on the chain are characterized by rectangular Young tableaux (of arbitrary height 
and width). The vacuum was proved to be of spin and formed by filled Fermi seas. 

In this paper, we carry on the study started in [31] with the calculation of the first excited states 
above the vacuum state. This construction generalizes the results obtained for the fundamental 
gl(N) spin chain [26,32,33]. The plan of the paper is the following. We first remind in section [2] the 
results obtained in [31]. Let us stress that the rest of the paper heavily relies on the results given 
there. In section [31 we consider hole excitations above the vacuum state, and the t corrections these 
excitations induce on the densities in the thermodynamical limit. We also compute the form of some 
hole excitations valid for any type of spin chains, as well as hole configurations corresponding to 
'small' excitations, and valid for a large class of spin chains (including all s/(2) spin chains and gl(N) 
alternating ones). More general excitations where new strings are added to the vacuum configuration 
are dealt in section HI The \ corrections to the densities are computed (theorem |4.ip and we give a 
form for the excitation leading to a state with 'small' (but non-trivial) spin. The Bethe equations 
relating the Bethe roots for holes and the ones for new strings are also computed (theorem 14.2ft . In 
section [5], we give a general form for a L -local Hamiltonian and compute the corrections to the 
energies for the excitations presented in sections [3] and HI Finally, appendix El is devoted to the proof 
of theorem 14.21 
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2 Notations and summary of previous results 

We give, in this section, a review of the previous results such that this article be self-contained. 



Bethe ansatz equations and string hypothesis 

In the papers [12,26,27], the Bethe ansatz equations have been established for spin chains where each 
site may carry a different representation of gl(N). To solve these equations, it is usual to use the 
string hypothesis which states that the solutions gather into Vm strings of length 2m + 1 (to G |Z + ) 
of the following form 

X^ k + ia, a — — to, — m + 1, . . . , to and j — l,...,N — l (2.1) 

where k = 1, . . . , Vm and \^ k , the center of the string, is real. Moreover, we will follow the lines 
given in [31] and consider representations which have a rectangular Young tableau. They are thus 
characterized by two integers (ae,ji) for the site £ = 1, . . . ,L, where an is the width of the Young 
tableau, and j e its height. The corresponding sl(N) spin is S — (0, . . . , 0, ae, 0, . . . , ). 

je-l N-l-j e 

Under this restriction, the BAE for the centers of strings can be rewritten as 

^^ ] ^%)-^Q%= (2-2) 

E E *-i m) (AS fe - ^) + e * { r\x m ] k - aS) + E **?\>i2* - 

v e\i+ \ <=i t=i i=i 

where the parameters Q^ k are quantum numbers. We have introduced the functions: 

$M (A) = ^ y? p+2a (A) , p G Z + , to G -Z + (2.3) 



m+n j 



$™(A) = ^ 2m+2 „+ 2 (A) + ^ 2 | m _ n |(A) +2 ^ ^ 2a (A) TO,nG-Z + (2.4) 

a=\m— n\+l 

m+n 

$5' n) (A) = - E VWA) (2.5) 

a=\m—n\ 

V? P (A) = 2Arctan , p G Z, p ^ and y9 (A)=0. (2.6) 

A state is completely characterized by the data of the numbers of strings, Vn \ and the quantum 
numbers, Q^ k - We can show that these quantum numbers are constrained by [31]: 

QZ e [-Q m,max i Qm,max ] (2-7) 
1 L 



= 1 



The valence, Pm\ is the number of allowed quantum numbers i.e. Pm = 2Qm,max + 1- 
Regular spin chains: definition and notations 

We want to consider the thermodynamical limit (i.e. L — > +00) but it seems impossible to deal with 
an infinite number of different representations. Thus, following the computations done in [31], we 
will work with L -regular closed spin chains, that is to say periodic spin chain (of length L) such 
that the representation (the gl(N) spin) at site £ is the same as the one at site £ + L , keeping L 
finite when L — > +00. Of course, for a L -regular spin chain, there is at most L different Young 
tableaux. However, among a sub-chain of length L , the same representation may appear several 
times. In this sub-chain, inequivalent representations may differ by distinct values of je and/or a^. 
We define C < L , the number of representations which have different values of d£, and called a a , 
a = 1, . . . , C these values. Accordingly, we define the ordered set: 

M = jn Q = Qq ~ 1 s.t. 1 < a < £ j C ^Z> with m < n 2 < ... < n c . (2.8) 
As a convention, we will also introduce no = — \ and nc+i = +00. The complementary set of M is 

(2.9) 



Jf= iz> \A/"= |m G iz> s.t. m^A/"| 



We also introduce the sets of indices defined by: 

I a = {£e [1, L ] s.t. a, = d a } , Va G [1, £] (2.10) 
such that ^ ^ 

E (■■■)< = E E (■■■>■ (2- 11 ) 

£=1 a=l £'e/ Q 

The cardinal corresponds to the multiplicity of a a within a subset of L sites. We define also 

which corresponds to the multipicity of the representation (a a ,j) within a subset of L sites. We 
thus have the property 

c N-l 

L ° = EE^'- ( 2 - 13 ) 
0=1 j=i 

We will also call J the greatest common divisor (gcd) of the J a /s. In most of the cases (i.e. as soon 
as a representation appears only once in the subchain of length L ), J is in fact equal to 1. 
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Bethe equations for the vacuum state 



Now, we can write down the Bethe equations in the thermodynamical limit for different states and, 
in particular, the vacuum state which is, by definition, the state such that Pm — Vm vanish for any 
m and j. In [31], we prove that this defines the state uniquely and one can show that it has spin 
zero (i.e. it is a trivial gl(N) representation). This corresponds to the following choice of parameters 



„0) 



NL, 



min(j, j e ) (n - max(j,ji)\ for 1 < a < £ and = for n G TV (2.14) 



°iei, 



Q% = k - + 1) for k = l,...,u$ andneAf . (2.15) 

We will consider the thermodynamical limit L — > oo, keeping Lq finite, of the Bethe equations 
(12.21) for a Lo-regular spin chain and for the vacuum state. Then, we obtain the following set of 
equations for the densities, a™ (A) (m G A/") of the center of the strings {A^\, | k — 1, . . . , oo}: 

/ dXat 1] (X)^T\^-X)+ / dA^)(A)^(A -A) 

„__., J-oo J -co 

+ r dXa^\X)^T\^-V) = -^^(^) + ^Y,(Y, 6 ^)^ ) ^o) (2.16) 
VA G ] - oo , oo[ , V? = 1, . . . , N - 1 , Vm G TV 
where ^(A) are derivative of $(A). 

Densities for the vacuum state 

To solve this set of equations, we perform a Fourier transform, with the following choice for the 
normalisation 

1 



/(p) = - J J*fWdX (2.17) 

and encompass them in a matrix. Finally, we obtain the following form of the BAE: 

27r§(p)£(» = A(p) (2.18) 

where the (N - 1)£ x (A - 1)£ matrix $(p) = -A(p) <g> and >l(p) is a (A - 1) x (TV - 1) 

tridiagonal matrix with the non- vanishing entries 

[A(p)] n = 2 cosh f M J and [A(p)] jj+1 = -1 = A(p)j+y (2.19) 

and \l/_i(p) is a £ x £ matrix such that [\I/_i(p)] a)i g = (p) . We have also introduced the 

(N — 1)£ vectors 

N-l c , N-l c 

£(P) = E E-S ef - 1} ® e(f) and A(p) = [EEE ^^(p) ^ ® # (2-20) 

j=l a=l j=l a ,p=llel0 



where ef is the canonical basis of C p . Then, inverting (12.181) and performing an inverse Fourier 
transform, we get the densities of the centers of strings for the vacuum state: 



sin 



ff W(A) = -±-Y y ~~^' NJ o Ka<£ (2.21) 

u eei a q=(\k-j t \+i)/2 cosh [j^^J ~ cos y-jj-J 

These densities allow us to compute physical quanties such as the vacuum energy. 

We will deal with sums that runs on integers or half-integers, with an increment which can be 
either 1 or |. To avoid confusion between these two types of sums, we will denote with a prime, 

, the sums that increment by | step, keeping the usual sum, f° r the ones with step 1. Hence, 
we have for instance the relations (for m integer): 



m-i 



£' = £+£ and £ = £+£ (2.22) 



n=0 n=0 n= | n= i n=\ n= i 



The aim of the present paper is to consider the first excitations above the vacuum state (as defined 
above). There will be of two different types: holes, which correspond to a configuration where some 
of the strings defining the vacuum state have been removed, and new strings which may be added 
once holes have been done. We will treat the BAEs for these excitations, that is to say gives (in the 
thermodynamical limit) the first corrections (in -jr) to the densities. We will also compute corrections 
to the energies and give 'local' Hamiltonians. 



3 Excited states: holes 

In this section, we consider the first excitations above the vacuum state. They are constructed from 
the vacuum by modifying the root distribution. For such a purpose, we will create holes (i.e. remove 
some strings) in the filled seas of the vacuum, and possibly add new strings. We first deal with holes 
created in the j th sea, suppressing some strings of length a p , p = 1,2,.... The parameters a p have 
to be some of the ai parameters defining the representations of the spin chain (i.e. ^ G M). 

3.1 Valences and spin 

Let {vn } be the configuration of the vacuum. The configuration corresponding to an excited state 
reads, for n £ M and 1 < j < N — 1, 

i* ) = ^-A*? ) (3-1) 
where the sets {/in \n £ A/"}, 1 < j < N — 1, label the seas where perturbations to the vacuum 
configuration are introduced. Let us stress that v£* as well as \x£) are a priori rational numbers 
(depending on L, the length of the chain) while must be integers (they correspond to the number 
of n-string in the sea j). In the following, we will use 

v<$ = 2$ - ^ - m £ N, j = 1, . . . , N - 1 (3.2) 
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with the convention fim = finP = 0. 

The corresponding valences are given by 

pjJ) = vji) + ^ min(2n + 1, 2m + 1) v ( £ . (3.3) 

They correspond to the number of allowed quantum numbers, while iffl denotes the number of 
parameters which are indeed used. Thus, for n G TV, 

pOO = p<J) - pf = min(2n + 1, 2m + 1) v$ (3.4) 

m£Af 

is the number of admissible quantum numbers \ which are not used. They provide the number 
(for each sea and each length of string) of free parameters attached to the excited state under 
consideration, i.e. the number of holes (formed by strings of length 2n+ 1) in the j th sea. Thus, they 
must be postive. Of course, some of the numbers T>$ may vanish. The spin of such a state is easily 
computed: 

Sj= £ (2m + 1) t$ = 2>g> . (3.5) 

rneAf 

It is equal to the number of holes in the j th sea of the longest string. 

In order to determine the first hole excitations of the model, one should look for the values of Vm 
such that Vm i T>n and Sj are positive integers and Z>„ and/or Sj have the smallest value. Such 
task seems impossible to solve on full general grounds, so that we will treat two slightly different 
problems: (i) Find excitations that exist whatever the type of chain (ii) Find 'small' excitations for 
a subclass of spin chain. The two following lemmas correspond to these two problems. 

Lemma 3.1 The hole excitations given by, for 1 < a < C, 



,(i) 



<Ja,j 



7f Yl 5 iM > ( 3 - 6 ) 



where J a j and J are given in and fallowings, exists for any L -regular spin chain. They have 

spin and hole numbers given by (for m G M and 1 < j < N — 1): 

V m = ■7f'%28 j)jt wm(2m+l,ai) and Sj = -j^2S jtje a e . (3.7) 
J £=i J e=i 

These values generalize the ones obtained in [15], for the spin s XXX model, and in [33], for the 
gl(N) spin chain built on fundamental representations (see also examples below). Moreover, since 
Vnl are integers, equations (13.41) and f)3.5p clearly show that the number of holes T>n and the spin 
Sj for this state are also integers. It leads, from the equations (13.21) . to the following values (for 
1 < a < £ and 1 < j < N - 1): 

— Y m ™(j,je)( N - max(j, j t )) (3.8) 
(Lo/J) 



'1q 



NJ 

^2 min(j, j e ) (N - max(j, . (3.9) 



u £ei a 
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We deduce that the corresponding excitation exists when the length of the chain, L, is such that 
(L — (LqI J)) is a multiple of NL . Such a length exists because J is also a divisor of L , see eq. 
f)2.13p . The factor 1/ J may look weird at first sight. However, it may be explained by remarking 
that the Bethe equations do not depend on the choice of the order of the representation*! on the 
chain. Therefore, by reordering the representations, we may transform a L -regular spin chain to a -y- 
regular spin chain with the same excitations. This may explain the presence of the denominator J in 
the above formulae. These excitations correspond to the highest representation in the decomposition 
of the tensor product of the representations entering in the ^-regular chain. The Young tableau 
corresponding to this spin is the juxtaposition of the Young tableaux (in decreasing order of ji) of 
all representations of the ^-regular chain (see third exemple below and Figured]). 

We now turn to the case (ii). The aim is to encompass a wide class of spin chain (i.e. a wide class 
of possible representations on the chain), and in particular for sl(2) spin chains, it will encompass all 
types of representation, provided^] Lq > 1. For gl(N) spin chains, the representations entering the 
Lo-regular chain will have rectangular Young tableaux of same height jo and arbitrary width at. 

Lemma 3.2 For Lo-regular spin chains with representations such that jt = jo, VI, and such that 
(NLq) is even, the following values define 'small 1 hole excitations: 

I 1—1 when a = L — 1, L — 3, . . . 

The hole numbers and spin are given by eqs (3.4) and $3.5\) . 



It leads, from equations (13. 2p . to the following values (for 1 < a < £ and 1 < j • < N — 1) 



Vnl = y min(j,jo)(N-max(j,j )) (3.11) 
»fi = 2L 2NL^ L ° min ^,jo){N - max(j, j )) . (3.12) 

We deduce that the corresponding excitation exists when the length of the chain, L, is such that 
{L ± ^y 1 ) is a multiple of iVXo. It is obviously possible only when NLq is even. This is in particular 
ensured for any sl(2) spin chains and also for alternating gl(N) spin chains. 

Remark 3.1 Apart from the excitations introduced in lemmas 13.11 and 13.21 there are many other 
possible excitations. In particular, any mutliple of these excitations is allowed: 

Ji^ = k^, VmeA/", j = l,...,N-l, keZ + (3.13) 

where /j™ correspond to excitations of lemmas 13.11 or 13. 2[ For instance, the spin 1 excitation of [34] 
corresponds to k — 2 for an excitation of lemma I3.2L 



2 On the other hand, the order is crucial in the computation of conserved quantities deduced from the transfer 
matrix: their explicit form may be completely different even though they have the same spectrum. In the same way, 
the order does matter for the calculation of Bethe vectors and correlation functions. 

3 Whcn Lq = 1, one deals with spin s sl(2) spin chain, whose excitations are treated in lemma HOI 



7 



Remark 3.2 For spin chains with a given length such that the vacuum exists (as determined in [31]), 
Vn^ are integers, so that excitations must have integer fjjp. We will say that the corresponding 
excitations lie in the vacuum sector. It is the point of view developed for example in [10, 15]. 

Examples 

1. We consider a gl(N) spin chain where all sites carry the fundamental representation N_. Since 
the vacuum state is built on real strings only, it is the only type of strings one can suppress. Then, 



and we get 



= V N JJ - ; = , n > (3.14) 



Sj = V® = 2H® - - , (3-15) 

in accordance with the results given above. 

The excitation carrying the sl(N) fundamental representation S — (1, 0, 0, . . . , 0) corresponds to 
the values given in lemma I3~T1 

A* W = ^- (3-16) 

The number of quasi-particles is T>n = #n,o<5j,i> as expected. This 'first' excitation lies in a spin 
chain with length such that L — 1 is a multiple of N. We remind that the vacuum is lying in a spin 
chain with L multiple of N. 

2. For a si (2) spin chain with a spin s representation on each site, the vacuum is built on s — | 
strings only. Thus, we get 

v s _± = n and v n = , Vn G -Z , n/s . (3-17) 

2 2 2 2 

The resulting excited state has spir0 f = 2s /x and a number V n — As /i S ns _i of spin | 'elementary 
particles' [10]. The hole excitation of lemma [37X1 corresponds to // = | (i.e. i> = 1), so that one can 
get a state with spin s. This spin s excitation is obtained when L is odd, while the vacuum has L even. 

3. If one considers an alternating spin chain, where the sites 2£+ 1 carry an sl(N) representation 
defined by (ax, ji) and the sites 2£ carry a representation (02,^2)- We suppose that a\ < a 2 , so that 
at = 2ri£ + 1. Hole excitations have configurations: 



1/^'^ for n^ni,ri2 (3-18] 



^ = for A; = 1,2 (3.19) 

The spin is 

= 01 (2a*£? - rfr 1} - + «2 - /i^ 1 ) - A*SL +1) ) (3.20) 



4 The factor ^ is implemented to make contact with the usual spin used for sl(2). 
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The excitations of lemma 13.11 are given by 



ij) j(N-jk) . . 

Kl = ~t J < Jk 



N 

Vnl = Jj 3k < J 



k = 1,2 



They lead to non-vanishing hole numbers 

VS = (2n 1 + l){5 hn +5 hn ) 

2?gJ = (2m + 1)^ + (2n 2 + 1)<^ 2 . 

The spin S = (Si, . . . , Sjv-i) is given by 

5j = for j 7^ ji,j2 and = for k = 1,2 

corresponding to a representation with Young tableau given in Figured] (drawn for j 1 > j 2 ) 



(3-21) 



(3.22) 
(3.23) 



(3.24) 




Figure 1: Young tableau of the hole excitation for an alternating spin chain (ax, ji)-(a 2 , j 2 ) 
If we suppose furthermore that j\ = j 2 = jo, one can consider excitations as in lemma [372} 



1)h j(N-j ) 



j < Jo 



k = 1,2 



(3.25) 



They lead to a unique non vanishing hole number T>n} = iV(n 2 — n{) Sjj and have a spin S = 
(Si, SW-i) with 

Sj = for j 7^ jo and Sj = N(n 2 — rix) ■ (3.26) 
3.2 Thermodynamical limit 



For m G Af, the quantum numbers for an excited state are 



1 1 V - 

Q% k = k--(^+l) + -J2^<k-k% d 



k = 1, 



(3.27) 



where we have introduced the sign function 

sign(j) 



1 J > o 
•1 j<0 



(3.28) 
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The numbers k~jj d are integers characterizing the holes. They must fulfill the constraint 

1 < *S, < ■ (3-29) 

As before, we take the thermodynamical limit of the Bethe ansatz equation and we obtain (see 
e.g. [15]): 



Y.\r dx ( s ^ i)(a) + s « +i) ( a )) *5 ,n) ( A ° - a ) + r dAs « )(A) ^ m,n)(A ° - a) 

n£Af U-oo J-oo 

= -2vr Bg?(Ao) + i E 5(A - Agj + j- E ( E ( 3 ' 3 °) 
VA G ] - oo , oo[ , Vj = 1, . . . , N - 1 , VmGjV 

where Sm (A) is the density of (2m + l)-strings in the sea j for the considered excited state. The 
parameter X^£ d is the image of k^ d /L in the limit L — > oo. 

Performing a Fourier transform, the Bethe ansatz equations become 

2tt$(p) S(p) = A(p) - - A(p) , (3.31) 

where \l/(p), A(p) are given in section [2] and 

AT-1 £ »n2 

A G°) = EEE^Of 1 '®^ (3-32) 

j=l a =i cJ=i 
AT-l £ 

2(P) = EE^SW^r 1 ^^- (3-33) 

j=l a =l 

By linearity, the solutions take the following form 

S(p) = E(p) + ±%) (3.34) 



JV-l £ 



where S(p) gathers the densities for the vacuum (I2.20p and i?(p) = E] E/ Pnliv) e j N ^ ® e cf' > 

j=l a =i 

corresponds to the densities of holes. It is solution of the following linear system 

2tt§(p) R(p) = -A(p) . (3.35) 
Inverting the matrix ^(p), we get the solution 

Rip) = ^ (^(p)" 1 ® (^-i(p))" 1 ) A(p) (3.36) 
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where the inverse of the matrix A(p) is given by (see e.g. [32]) 

\P\\ ,^ f „•>!' 



sinh (^(N — max(i, sinh ^min(i, j) - ^ - ) 

sinh (ivM) sinh (M) 



C4(p)~% = " VTT, — • ( 3 - 37 ) 



The inverse of the £ x C matrix \&_i(p) is provided by the following lemma (we remind that we 
ordered n\ < n 2 < ... < n c ). 

Lemma 3.3 The inverse of ^/^i(p) is a symmetric tridiagonal matrix. Its diagonal elements are 
given by, for 1 < a < C 



OMp))" 1 = - sinh (M) sinh Q| (n a+ i - Wq-Q) 
V /aa sinh (l»l — n„,„i^ sinh (l»l — 



where, by convention, we setn = —\ andnc+i = +oo. We also used the limit s™h(^-bj — exp(6 — a). 
The non-vanishing non-diagonal elements read, for 1 < at < C — 1, 



($-*))"' - . hn s l" h( " ) = ($-*))"" 

V / a,a+l smh ( » \ (n a+ i — n a ) I V / a+l,a 



*+i sinh (n a+ i - n Q )) 
In £/ie particular case £ = 1, this matrix reduces to a number: 



(3.39) 



\-i Id p(bl(™i+f)) 

*-'<» - 3i " h (T) Bm :n B ir,, + 



2 ' sinh (m + i)) 

Proof: Direct computation on the product ^\&_i(p)j ^_i(p). 
Let us introduce useful functions, for 1 < a < £, 

27r „^^rN 



(3.40) 



i sinh f fiV — maxfi fc)) — 1 sinh ( mm(j, k)— 1 
ASM = - i U , ,., 2 , J 1 2 ' (3.41) 



sinh (^y) sinh ((n a+1 - n a )\p\^j 
and 

?s (?) = (?) — v r (3 - 42) 

smh \ {n a - n a _i) |p|J 
where the convention no = —1/2 and n^+i = oo has been applied. 

Theorem 3.4 Let us consider a L -regular spin chain based on gl(N). For the excited states corre- 
sponding to holes in the filled seas, the densities are, for 1 < j < N — 1 and 1 < a < C, 

*2(A)=aJ2(A) + ipg)(A) (3.43) 
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where <7n2(A) are the densities \2.21\) for the vacuum and the corrections at order 1/L are 



N-l 



fL a — l Yi-ct rt a-f-l 



m = e 



k=l 



(k) 



d=l d=l d=l 



(3.44) 



T/ie functions h&(\) and g^{\) are the inverse Fourier transform of ( 3.41 ) and ( 3.^2^ respectively. 



3 

By convention, we used VnJ = *D_h 2 = and T>n} +1 = = 0. 

The corrections to the densities, for the particular case C = 1, reduce to 

N-l V\{ 

pS(a) = E ( 3 - 45 ) 

k=l d=l 

where gjk{\) = gf k (A) (for C = 1, n = -1/2, n 2 = oo). 

Proof: We give the proof only for a = 1. The other density corrections are obtained similarly 
Projecting the linear system (I3.36f) on its first component and using the lemma 13731 we get 

N-l /< < \ 

pt? (P) = E E tk (P) exp(ipA^J - E (?) °P(^ ( 3 - 46 ) 

k=l \d=l d=l J 

Performing an inverse Fourier transform and using the conventions n$ = — ~ and = 0, we get 
the result fl3~44"|) for a = 1. ■ 

Let us remark that the explicit expressions of h^ k (X) and gf k (X) are not needed, since physical 
quantities are computed via their Fourier transform (see for instance section [572]) . 

3.3 Examples 

1. For a gl(N) fundamental spin chain, the corrections to the real root densities are 
1 sinh((JV-maxO',A))M) sinh(minO,*)M) 

#'<p) = - s E E — * -. w 2 ' A iJ - ^ + " x - ^ 

i-i -i-i sinh ( sinh (y 



where we noted X d for Aq d . 



For a pZ(2) spin s chain, the corrections to the root densities of length 2s are 

\P\ 



where A^ = ,. 



P.-i/a(p) = -o- E • ui n J I h eS ' Pl+tPAd ( 3 ' 48 ) 
27r ^— ^ smh(|p|J smh(s|p|) 
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3. For a gl(N) alternating spin chain with ni < we get 



and 



N-l 



( \v\\ ( \v\ 

sinh \ {N — max(j, k))~^~J s i nn y m i n (i) k) 



m = - 1 E 



2vr ^ . , f nT \p\ 



k=l 



sinh (^N-^J sinh ^(n 2 — ni) \p\ 



N-l 



sinh((n 2 + |)b|) - tp ~ x ^ a _ J* ^ . 
sinh((n 1 + i)|p|) ^ ^ r 



sinh y(N — max(j, k))~^~J s i nn (^ m i n C/\ ^)~2/ 



2vr ^ . . / A Jp| 



sinh ^N— ^ sinh ^(n2 — rii) |p| 
— ^ e' pA "i' d + exp((n 2 — e ipA ™2> d V (3.50) 



4 General excited states 

We now turn to a more general excited state. It corresponds to the case where holes are created in 
the filled seas and new strings added in other seas. 



4.1 Valences and spins 

The configuration of an excited state is characterized by 

= z/<f) - jig) f or n E M and u { n j) > for n E 77, (4.1) 

where we kept the notation Un , n E Af, for the vacuum configuration. The set 77 (complementary 
to M) has been introduced in (12.91) . Obviously, if iffl = 0, Vn E 77, we recover the previous case 
with holes only. 

The corresponding valences and spins are given by, for n E 

PP = tiP + E' min ( 2n + !> 2m + 1) (2mS? - &- x) - ^ +1) ) 

me A/" 

- min ( 2 ™ + 1, 2m + 1) (2*g> - Jg" 1 ) - ^ +1) ) (4.2) 

me AT 

5, = £'(2m + 1) (2/ig) - /4T 1 ) - /^ +1 >) - £'(2m + 1) (2Zg> - ^ - • (4-3) 

Obviously, demanding the valences and the spin to be positive imposes contraints between the number 
of holes, V { i } = p£> — Vm [m E M), and the number of new strings, v£ (n e77). We study in more 
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details particular excitations in section 14.41 

The new strings of length 2n + 1 with n G A/* are characterized by quantum numbers {Q n k | 1 < 
k < Vn } and by the Bethe roots {X^ k | 1 < k < Vn (these sets may be empty). These quantum 
numbers are bounded, namely, for n G AT, 



1 - / ^ P {j) - 1 

" <^,i<^V^- ( 4 - 4 ) 



Let us emphasize the difference between \ (n G AT), Bethe roots for the new strings (the ones 
added to the vacuum configuration), and X^' k (n G A/"), missing Bethe roots in the filled seas (holes 
in the vacuum configuration). 

4.2 Densities for a general excited state 

In this section, we compute the densities for the general excited states in the thermodynamical limit. 
This is done in two steps: first, the Bethe equations for n G M provide the correction to the densities 
in terms of the new Bethe roots \^ k , n G AT (see theorem 14.11) . Second, the Bethe equations for 
n G A/* couple these new Bethe roots with the holes in the filled seas (see theorem 14. 2p . These 
relations depend on the choice of the quantum numbers {Q^ k \ n G A/", 1 < k < Vn^}- 

Using again linearity of the Bethe equations, the densities for these excited states rewrite: 

t?(A) = 4 J) (A) + ^^(A), n G A/" , (4.5) 



'(A) = Z E^ A - A ^)' ( 4 - 6 ) 

where Sn (A) = er^(A) + £/>n (A) are the densities when only holes are created in the vacuum 
configuration, as given in theorem 13.41 while Cn ^(A) are the corrections due to the new added strings. 
We have assumed that the number of these new strings remains finite even in the limit L — > oo, 
hence the form (14. 6 1) . 

From the Bethe equations for m G M and performing a calculation similar to the one of section 
|3l one shows that Cn\\) must satisfy the following equations 

e' { r dx ( c ^ i)(a) + c - +i)(A) ) ¥t){x ° - a) + r dx c - )(A) ^ m,n)(A ° - a) ) + 2n c - (Ao) 

w — oo J —oo ) 



n£Af 



oo 



£' E *5 ,n) (A - xgr 1 ^ + E *2 m,n) ( A ° - *$) + E *5 ,n) (Ao - a^ i} ) = o (4.7) 

neJT ( <=i '=1 «=l J 

VA G ] - 00 , 00 [ , Vj = 1, . . . , N - 1 , Vm G AT 

Before solving these equations, we introduce the maps 7 ± : 

I 7~(m) = n a G A/"U {n = -\,n c+ x = 00} 
m G N 1— >■ < ~ with n a < m < n a+ i . (4.8) 

I 7+(m) = n a+1 GjVu{fi = -§,nc+i = 00} 



14 



In words, the functions 7^ associate to any element in TV", the two closest numbers belonging to 
7VU {-i,oo}. 

We also need the following functions, for m G 77, 



sinh (p \m — 7 ± (m)|) 



C m(p) = ■ 7-7 ) + , ; , ^ • (4-9) 
smh [p (7 + (m) — 7 (mJ)J 



Their inverse Fourier transform is 

71 \m — 7 ± (m) 



sin 



7+(m)-7 _ (m) , / vrA \ / 7T m-7 ± m 1 x 

' v ; ' v 1 cosh ——- -— + cos ' 



7+(m) — 7 (m) / \7 + (m) — 7 (m) 

Let us remark that, for m G A/* finite and such that 7 + (m) = 00, these functions reduce to 

^( A ) = T2X7 -( \v and <( A )=° ( 4 - n ) 

A z + (m — 7 (m))" 1 

Now, we can give the corrections to the densities for a general excited state. 
Theorem 4.1 The densities for a general excited state are given by ( |^.6p and by 

<ff(A)=*?(A) + -]rcg>(A), nGAA (4.12) 



where Sn (A) are gruen m theorem 3.J^, The corrections Cn (A) are gruen /or 1 < a < £ and 
1 < J < iV - 1, 



c S(A) = -^ E + E E^( A - A S) ( 4 - 13 ) 

\ n a -\<m<n a 1=1 n a <m<n a j r \ £=1 J 

We remind that uq = —1/2 and nc+\ = oo. By convention, we also set u)q(X) = 0. 

Proof: Performing a Fourier transform of equation (14.71) and inverting the linear system, we get, for 
1 < a < C, 

1 C 

^2(p) = ~EE ex p MS) E (*-i(p)) a /-T m) (p) • (4.i4) 

71 mGA/" <=1 (9=1 



From the lemma I3.3[ the sum on (3 in the previous equation can be computed and we get finally 

ria-h ?iU) n a+ i-i ~U) 

= E' E s ™(p)expM£,)-^ E' E^W ex p^ A S)> ( 4 - 15 ) 

The final result is obtained through inverse Fourier transform of this last equality. ■ 
In theorem 14.11 corrections to densities have been expressed in terms of the parameters A 



(m G AT). These parameters are implicit functions of \^\, as stated in the following theorem. 
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m,i 



Theorem 4.2 The relations between the Bethe roots corresponding to holes in the filled seas and the 
ones corresponding to new strings are given by, for m G AT", 



V (J)_ 

7 (™) 



7"r (rn) 



7+M ,J (4-16) 

d=l <2=1 
r= 7 -(m)+i I <=1 <=1 «=1 



where the functions are a primitive of uo 



^m( A ) = 2 arctan 



tan 



it m — 7 (m) 



tanh 



2 7+(m) — 7~(m) / \2 7+(m) — 7~(m) /. 



7T 



A 



(4.17) 



and the functions Fa T ' m \\) are 

Fi r r } (A) = 



F^ m \\) 



r«(A) 



r+m— 27 (m)— i 

£ r S r) ( A ) 

g=|r— m| + | 

2m-l-27~(m) 

( A )+2 £ rW(A) 



(4.18) 



2m— 2 7 — (m) 



if m = r 



q=l 



r 



(r) 



r+m— 1— 27 (m) 



,+„,-+. „»)(*) + rf r l m| (A) + 2 £ r ^( A ) ? ™ ^ r 

q=\r—m\+l 

n q — 7 + (r) + 7~(r) ' 



-2 arctan 



tanh 



7T 



A 



tan 



2 7 + (r) — 7~(r)/ V2 7 + (?~) — 7 _ (^) 



(4.19) 



(4.20) 



Proof: Postponed in appendix El 



Let us remark that if m, r are such that m, r > ri£, the functions fi+(A), rjf^A) and 

F a (r,m) (A) (a = -1,2) reduce to 



fi+(A) = 2 arctan f 



A 



m — nc 



and 



, rJ m )(A) =2 arctan (-) = y? 2(? (A) 



(4.21) 



'(A) . (4.22) 

Notice also that F^ m \\) = F^ m ' r \\) when 7 (m) =7 (r), which is always the case in the relations 
where these functions are used. 



16 



4.3 Examples 

1. For gl(N) fundamental spin chain, we get AT = {|, 1, §, . . . } and j~(m) = 0, 7 + (m) = oo for any 
m £ A/. Thus, the corrections for the real Bethe roots reduce to 



oo vi i] 



c 



«(A) = - E I— |- (4.23) 



The constraints (14. 16f) become, for m G AT, 

v / I 



= £ *>-(A« "f)-E E *?r""" 3, (A« - AOr") (4.24) 

d=l r= i [ £=1 



We remind that stand for A^. 



For the particular = for any m > 1 (i.e. only strings of length 2 are added), we get 

rrt?) rrU) r:(i) 

^1 v \ v \ 

2,Qf t = E^(Aa-C'>)-E^(Aa-A^) + t^(A£-A««') 



f=l 



<i=l 



We recover a previous result computed in [33] . 

2. For a g/(2) spin s spin chain, we get AT = {0, ~, . . . , s— §, s— 1, s, s+|, . . . }. Thus 7 + (m) = s— ~, 
7~(m) = — | for < m < s — 1 and 7 + (m) = oo, 7~(m) = s — ^ for s < m. The corrections are 

, ,U) = fivi ""'^ ^'Y- 1 m ~ S + ^ (126) 

i=S tT 2s cosh(^^) + cos(x^±l) t,tl*( X - A™.') 2 + ("> - « + ^ 

The constraints become, for m = 0, |, . . . , s — 1, 
27rg m , fc = 2 53 arctan tan (^j) tanh [ E^^g) j _ £ F 2 (r ' m) (A m , fc - A P ,<) (4.27) 

d=l |_ ^ ' •■• ■ ' . l /=i 

and, for m = s, s + ~, . . . , 

2vrg m , fe = 2 53 arctan ( ^ ~ ^ ) - J^E $ ?~ S * m ~ S) (A m , fc - A r/ ) (4.28) 



__i <=i 

2 £^fe 



z> 1 



d=l \ 2 / r =s «=1 
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keeping the notation = A s _i d for the Bethe roots of the holes. 

3. For an alternating spin chain with ri\ < n 2 , we get 7~(m) = — ~, 7 + (m) = n\ for < m < 
ri\ — |, 7~(m) = ni, 7 + (m) = n 2 for m + | < m < n 2 — \ and 7 _ (m) = n 2 , 7 + (m) = oo for 
n 2 + \ < m. Thus, the corrections to the densities are 



2 sin f 



i *. i Bill 7T 

* m=0 <=1 + 2 cogh / ^^V, \ , cos 



n2 ~2 Sin (Tr ^+^-ni 

1 x — v x — > 1 \ n,2— ni 



9 ^ ^ r>„ - 



m=ni+ i £=i ™ 2 ni cosh f vr^^l + cos f tt "^" 2 -" 1 



and 



sin I 7r m+n2 - ni 



=ni+i COsh 7T ^ + COS 7T- 

1 ri2— ni / V n-2,—n\ 



TV 



m — n 2 



- m ^ +| tr(A-AS) 2 + (m-n 2 p 

The constraints between the roots are still given by the theorem 14.21 Since it is a rather cumbersome 
equation, we refrain from giving its explicit form. 

4.4 Excitations within the vacuum sector 

We consider an excitation in the vacuum sector, as defined in remark I3.2t its hole configuration 
is given by {Jim , m G jV, j — 1, ...,N — 1}, where Jim = N /j, m ' and {^m} describes the hole 
configuration as given in lemma [37X1 To this configuration, if we suppose that rii > 0, one can add 
new strings of the following type. 
We first introduce (for m G TV 7 ): 



8 m ,n a ^nl where n a = max(m G J7 , m < n a ) . (4-31) 



a=l 



We remind that v$ = 2J&- - (m G 77) and u#> = 2 ^ - 



The set {/i™ , z/ir , m G jV , n <Ej\F , j = 1, . . . , N — 1} define a configuration with a spin 

= ^2(n Q -n a )^2 (4-32) 
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and valences 



P m = v m if vn<n x 

fa c 

a=/3i+l a=/3 2 +l 
£ 

a=l 



(4.33) 



Above, /?i and ^ depend on m and are such that 

< m < ^jSi+l an d ^/3 a -i < 171 — ™P2 (4.34) 
This construction generalizes the construction given in [10] for the spin s XXX spin chain. 

Example 

For the gl(2) spin s spin chain, the lemma I3TT1 hole excitation is given by v s _i_ = ^y^, leading to 
w n = 2S s -i >n , i.e. v s -i = 1- Then we get an excited state with spin ~ = 1, corresponding to the 
adjoint representation of sl(2). As already mentioned, we recover the results given in [10]. 



5 Hamiltonian and energies 
5.1 Lo-local Hamiltonians 

Definition 5.1 An Hamiltonian is called Lo-local when it possesses only interactions between neigh- 
bours separated by at most L — 1 sites. 

This definition is quite natural. It is used in the following: 

Conjecture 5.2 For any L$-regular spin chain there always exists an Hamiltonian which is Lo-local. 

This conjecture is supported by the theorem 15.31 given below. To formulate it, we first introduce 
local representations of gl(N). On each site I we denote the representation of the gl(N) algebra by 

TV 

£14 = ^ (S e$ (5.1) 

where the indices label the auxiliary and/or quantum spaces, the generator acting as identity operator 
in the other spaces. Accordingly, the Yangian representation at site £ will be Ti^(A) = A Ijv<g>I+z£i^. 



19 



Theorem 5.3 For a gl(N) spin chain, L^-regular, with a fundamental representation at sites 1+gLo, 
q G Z +; the Hamiltonian H = — it(O) -1 1'(0) is L -local and takes the form: 



L> 



e it 



q=l 



1=2 



l+qL ,Z+{q-l)L 



1 / —> 



~ s n Si +i L ^ 

£ 1= 2 \j=2 

where the total number of sites is L = pLo . 



Y[Sl + ( q -l)L ,£+(q-l)L ] Pl+(q-l)L ,l+qL 
\ 



e=2 



^1 



(5.2) 



/ 



Proof: We consider a gl(N) spin chain, Lo-regular, with a fundamental representation at sites 1 + qLo, 
q G Z + . We define the Hamiltonian as 



= -it(0)~V(0) with t'(X) = -^-t(X) 



(5.3) 



Starting with t(X) = tr (T i(A) . . .Tql(X)), we decompose t'(0) into three terms, depending whether 
the derivation applies on T i(A), on T ^(X), 1 < £ < L + 1 or on T ^(A), L < £. After multiplica- 
tion by i(0) -1 , the first two terms can be easily computed, while in the last one, we recognize the 
Hamiltonian of the same type of spin chain, but where the first L sites have been suppressed. Then, 
we get the recursion formula for the Hamiltonian: 



H < L +1,L +2,...,L> + ( ]^[£l+Lo/ 

t=2 



<1,2,...,L> 



1+L ,£ 

e 1= 2 \e=2 

This recursion is easily solved, and we get (15.21) . 



IT 

1=2 
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Pi 



1+L 




1 / 



e=2 



(5.4) 



Of course, depending on the representations sitting at sites £ ^ 1 + qL , further simplifications 
may arise in the expression (15.21) . see examples below. 

Note that once the representations of the spin chain are fixed, £\ i becomes a matrix, which can 
be easily inverted. In a more algebraic way, one may compute (£i^) _1 using the quantum comatrix 
in the Yangian Y(N). 

When no fundamental representation occurs on the chain, one has to deal with a fusion procedure, 
in order to get auxiliary space with higher dimension (see e.g. [28,35]). It is reasonable to think that 
to get a Lo-local Hamiltonian, the auxiliary space has to be isomorphic to one of the representation 
on the chain. Then, the transfer matrix will be of different type. However, the Bethe equations (and 
thus the results presented above) will remain unchanged, since all fused transfer matrices commute 
among themselves. 
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Examples 

1. For gl(N) fundamental spin chain, one recovers the usual spin chain Hamiltonian 

L 



H = Y,P^+i (5-5) 



e=i 

2. The spin-s gl(2) chain does not enter the hypothesis. However, the existence of a local (1-local) 
Hamiltonian is known [15,36]. 

3. For an alternating spin chain (with L = 2p), with fundamental representations on odd sites, one 
gets an Hamiltonian with next-to- nearest neighbours interactions (i.e. 2-local in our terminology): 

L/2 

H = {(^2£+l,2i) 1 ^-1,2^^2^-1,2^+1 — ^2e+l,2l} (5-6) 

If furthermore, one takes a gl{2) spin chain, the matrices 8 and £ _1 are easy to compute. Indeed, 
in a sl{2) representation of spin s (square matrices of size 2s + 1), one has 



n=l n=l 

2s 2s+l 

= Vn(2s + 1 - n) E^ n ; tt s (1) = £ E$ = I 2s+1 (5.7) 



n=l n=l 



where E\f are the (2s + 1) x (2s + 1) elementary matrices with 1 at position and elsewhere. 
In this representation, one has 

[e 3 ,e±] = ±e ± ; [e + ,e_]=2e 3 and tr{e+ e_) = tr(e 2 3 ) = s ( g + 1 K 2s+1 ) ( 58) 

Hence, we get 

2 £2^+1,2^ — h ® Wi + 4 (<r 3 <g> 7r s (e 3 ) + - <g> 7r s (e_) + (j_ <g> 7r s (e+))) = I + 4X S (5.9) 
where we have introduced the Pauli matrices 

*s=(g A) ; ^+=(U J) ; 0) ■ ( 5 - 10 ) 

Due to the properties of Pauli matrices, we have (whatever the value of s) 

X s 2 = il 2 ®C s -ix s (5.11) 
where C s is the value of the Casimir operator in the spin s representation: 

C 8 = 7r s (e£ + ^(e_e + + e+e_)) = s(s + 1) I 2s+1 . (5.12) 
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This implies that 

/ 1\ 2 1 
(£2e+i,2t) 2 = [s + -J I 2 ®Wi ($u+i&)~ = — 2^+1,2^ (5.13) 

\ S + 2) 

Using the notation 

aW- a ■ S 2i = ® vrf )(e 3 ) + <xf ^ ® Trf >(e_) + a^^ ® Trf >( e+ ) (5.14) 

where cr^ ^ and 7r« (e3,±) are representation of at site 2i — 1 (spin |) and 2£ (spin s) resp., 
and obvious generalisations of this notation to other cases, we can rewrite the Hamiltonian as 

-1 L/ \ 

H = — —2 V, \ s ( s + !) + 4s ( s + 1) ■ oWi - • Sa£ - oW-i • <?2£+i 

2(s+2j <=i 

—4 (S^ ' ? 2£+l)(o ? 2£-l - S 2 £) ~ 4 (5*2£ ■ <T2l+l) {&21-1 ' ^Wl) — 4 {3 2 t-l ' SW) (<?W-1 ' ^Wl) 

-16 (S 2 * • a 2 ^+i)(oWi • Su){vu-i ■ ^Wi)} (5.15) 

which describes an alternating gl{2) spin chains with spins \ and s and next-to- nearest neighbours 
interaction. 

5.2 Energies 

Motivated by the result given above, we focus here on the case of Lo-regular spin chain containing at 
least one fundamental representation. This implies that n\ — (we remind that the n^s are ordered). 
In this case, the energy (eigenvalue of the Hamiltonian ( 15.21) ) per site for a state characterized by 
the densities ti j) (A), j = 1, . . . N - 1, n E ~Z + ( as given in ( 14. 5ft and ( 14.61) ) reads 

n „ x -v /• +0 ° „ m/ , x 2m + l 

E~E + = I ^(A) x2 + {2 _ m±1)2 (5-16) 



meZ+/2 

where E + is the energy of the pseudo-vacuum 

c 



^ + = -l-EE^— ( 5 - 17 ) 



In particular, for the vacuum state (i.e. when densities are given by (I2.2ip ). one gets [31] 

e - e + = — y y L{ *°- jt+1 + 1) - $$*±l±) } (5.i8) 



where ^(x) = ^lnT(a;) is the Euler digamma function. For the excited states studied in this paper, 
the correction to this energy are given by the following theorem. 
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Theorem 5.4 For a Lo-regular spin chain with, at least, one fundamental representation, the cor- 
rection at order 1/L to the energy of the vacuum, due to the general excited states introduced in 
section^ is 



AE= 2 4YY Sm W (5.19) 

Proof: The corrections at order 1/L to the energy, due to the modifications of the densities in the 
seas m G Af, are, respectively, 



A^=E / dApg?(A) and A S £=]T / rfA c«(A) 



2m + 1 



A 2 + (^±!) 2 

(5.20) 

which can be computed via the Plancherel's theorem 

, r+oo , p+oo 

A h E = 2nJ2 / dppfflfye-^W and A S E = 2ir £ / dp?$(p) e -^ |p| . (5.21) 
The corrections to the energy provided by the densities (14. 6p of new strings read 

*™ = £' E Try— — — — — • (5.22) 

meA 7 «=1 ^ A m,^ I 2 J 

Using the explicit forms (13.441) of the functions p$(p), we can show that, in the expression of A h E, 
the coefficent of exp(ip\^ d ) identically vanishes for all n a but n\. Thus, one obtains (reminding 
that rii = 0): 

V - |P " - S inh((Af - fc)M) ,„™ 



fc= i d=i ■/-«> smh(^) 
Similarly, using the explicit forms ( 14. 13ft of c£v (p), we can simplify A S E: 

A S E = d P e"^ |p| e ipA ^ (5.24) 

_-T7 «_1 J —OO 



Then, performing an inverse Fourier transform of (I5.23P and ( 15.24ft . we get the results (15.191) by 
remarking that A S E = —A n E. m 



We would like to emphasize that although the studied spin chains are quite general, the final 
expression (15 . 1 9[) for the corrections AE to the energy are very simple. Indeed, these corrections 
have the same shape for any Aq^, and can rewritten as: 



A^EE^'S) where jM(A) = | , fc7r ■ (5.25) 
t^i^i N cos (lv ) - cosh ( lv A ) 
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Therefore, the holes in the filled seas of real strings can be interpreted as particle-like excitations 
with energy 5^(X ( k) d ) and rapidity \™. 

Moreover, since only real holes appear in the expression of AE, one can say that, at order i the 
holes created in filled seas of non-real strings have a zero energy. 



A Proof of theorem 4.2 



The parameters X^\ satisfy the following constraints which are provided by the Bethe equations: 



Y!\ d\ (fg-D(A) + tg +1 >(A)) SLT ^ - ^ + / dX #>(A) <S>t' n \*i] k ~ A) 

~0) ~0+i) 

*=i 

+rE(E«**)* ! ? ) < J & ) -T Q & 

a=l ^eJ„ 



(A.l) 



where m G A", 1 < fc < i/m and 1 < j < N — 1. To simplify this relation, we replace the densities 
t£°(A) in the l.h.s. of (jATTp by their values obtained from (g33D and (02]) : 

t? (A) = 4° (A) + ±p?(A) + ^cSP(A) . (A.2) 
The three terms of the above sum are reduced in the three following lemmas: 
Lemma A.l For me AT and Xq G R, we have the identity 

£' { r dX (aS^HX) + o-^(X)) $5' n) (A - A) + r dA<#>(A) $f-)(A - A)| = 

rE(E^O^ )(Ao) (A - 3) 

Proof: We derivate w.r.t. Ao the l.h.s. of (1A.3jl . perform a Fourier transform and use the explicit 
form of the Fourier transform of the vacuum densities 02.21 j) . Then, we remark that it is equal to 
the Fourier transform of the derivative of the r.h.s. of ( 1A.3I) using, in particular, 



hj-i t k(x) + hj + i t k(x) — 2 cosh(x) hj } j.{x) = —Sj^ sinh(A:r) sinh(x) 
where hj^{x) = sinh (x (N — max[j, k])) sinh (x min[j, k]) . (A. 4) 



This proves the equality up to a constant, which is fixed by considering the value Ao = in the 

3 



equation, and remarking that the densities are even functions while $^ m ' n ^(A) is odd. 
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Lemma A. 2 For m G A/* and for Aq € R, we /iai>e i/ie equality 



E' I r rfA +^' +1) (A)) $5' n) (Ao - A) + r dApW(A) $r-)(A - A) 



n£Af 

7 ( 



7~r (m) 



7+(m),d' 



V 



d=l 



d=l 



(A.5) 



Proof: Performing a derivation w.r.t. Ao and a Fourier transform of the l.h.s. of (1A.5j) . we get, using 
expression (I3.44p . the equality (1A.5[) up to a constant. Considering the limit A — > oo, one shows 
that the constant vanishes. ■ 



Lemma A. 3 For m 6 AT, an<i /or A G K, we have 

e' i r ^ ( c ^ i} ( a ) + c n +i) ( A )) $ 5 ,n) ( A ° - a ) + r rfA c « )(a) ^ m,n)(A ° - a) ) = (a - 6) 



v / I 



-EE $ -i m) (Ao - A^ 1} ) + E $ ?* m) ( A o - A $) + E $ - 9 i m) (Ao - A^ 

geA/ 7 [ «=l <=i <=l 

+ E' E ^(Ao - A^ 1} ) + E^' m) (Ao - + E ^(Ao - A& +1 >) 

r=7 -( m )+i [ *=i *=i ^=1 

where F_{ m \\) and F^ r (X) are given by an d 



Proof: Performing a derivation w.r.t. Ao and a Fourier transform of the l.h.s. of (1A.6j) . we get, using 
expression (I4.13P 



I 



M- 1 ) t;U) 

i^g iSq 



E E ^(p) exp(^A^ 1} ) +E § ?' m) (p) exp(zpAg) + E ^T\v) exp(ipA< 
qeir { e=i <=i <=i 

E' E ffi m) w ex P (zpA^ 1} ) + E ff ,m) (p) ex p(^S) + E £i"°(p) «p(^ +1) ) 

r= 7 -( m )+| [ <=1 i=l e=l 



where 



c(r,m)/ \ _ sinh(|p|(7+(m) - max(r, m))) sinh(|p|(7 (m) - min(r,ro))) 
/ - 1 W sinh(|p|/2)sinh(|p|(7+(m) - 7"(m))) 1 ' J 

_p> m )(p) = ~ 2sinll (bl(7 + (^) -max(r,m)))sinh(|p|(7~(m) -min(r,m))) ^ 

tanh(|p|/2) sinh(|p|(7+(m) — 7~(m))) 
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To obtain the explicit forms (14.181) and (14.191) . one can perform the inverse Fourier transform by the 
theorem of residue and the integration. Remark that it is simpler to check these forms through a 
derivation and a direct Fourier transform. This fixes the equality up to a constant. Considering the 
limit A —>■ oo, one shows that the constant vanishes. ■ 



Finally, using relations (14.51) . (IA.3I) . (IA.5j) and (IA.6j) to simplify Bethe equations (IA.1I) . we obtain 



the relation ( 14.161) . which ends the proof of theorem 14. 21 
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